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ABSTHACT 


The  computation  of  parameters  for  symmetric  laminates  in¬ 
volves  series  of  tedious  calculations  for  both  in-plane  and 
flexural  properties.  In  order  to  alleviate  the  magnitude 
and  conplexi*y  of  calculations  for  coeplicated  symmetric  la¬ 
minates,  the  program  of  this  report  mas  designed  for  a  com¬ 
mercially  available  microcomputer,  the  Hewlett-Packard 
HP-87,  and  is  based  for  the  most  part  on  program  logic  de¬ 
veloped  by  the  Air  Force  Bright  Aeronautical  Laboratories 
for  the  TI-59  handheld  programmable  calculator.  Program  lo¬ 
gic  is  explained  in  detail  with  flowcharts  and  a  full  list¬ 
ing  of  the  program  is  included.  A  description  of  the  pro¬ 
gram  logic  provides  the  user  with  a  comprehensive 
explanation  of  program  operation. 
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Chapter  I 
INTRODUCTION 

The  program  of  this  report  is  based  an  class  work  per- 
foraed  during  the  suaaer  Quarter  1982  at  the  Naval  Postgrad¬ 
uate  School  and  taught  by  Dr.  Stephen  R.  Tsai  of  the  Air 
Porce  Wright  Aeronautical  Laboratories.  The  original  pro- 
graa  was  based  on  reference  1  and  was  written  in  order  to 
coapute  and  observe  a  variety  of  in-plane  and  flexural  par¬ 
ameters  calculated  for  syaaetric  laminates.  Upon  completion 
of  than  program,  however,  it  was  enhanced  by  the  addition  of 
subroutines  and  displays  to  compute  the  max  stress  and 
strain  strength  ratios  in  addition  to  quadratic  strength 
ratios. 

Chapter  II  explains  program  logic  based  on  the  flow 
charts  of  Appendix  A.  The  reader  will  find  that  no  special 
programming  techniques  were  required  and  all  calculations 
are  straight  forward. 

There  are  two  basic  advantages  to  having  a  program  such 
as  this  on  a  microcompter.  First,  although  a  microcomputer 
such  as  the  HP-87  is  obviously  not  as  portable  as  a  hand¬ 
held  calculator,  the  user  will  note  distinct  advantages  to 
being  able  to  view  labeled  output  and  matrices  several  lines 
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at  a  tiae  as  opposed  to  single  line,  unlabeled  output.  Sec¬ 
ondly,  the  prograa  code  is  very  easy  to  read  and  aodify  if 
necessary. 

This  prograa  is  not  intended  to  be  all  encoapassing,  but 
it  should  be  of  assistance  to  users  designing  syaaetric  la- 
ainates.  Osers  are  encouraged  to  aodify  the  prograa  to  suit 
their  own  needs  and  to  publish  addendua  and  erratta  to  this 
report.  Corrections  and  nodifi cations  should  be  addressed 
to  Dr.  Tsai  at  APWAL/HLBB,  irigh t-Patterson  AF3,  Ohio. 


Chapter  II 
PBOGBAH  LOGIC 

Even  though  this  prograe  has  been  written  fcr  the  HP-87, 
modifications  to  fit  this  program  on  other  aicrocoeputers 
should  not  be  difficult.  The  only  really  unique  aspect  of 
the  progras  occurs  during  the  initial  parts  of  the  prograe. 
Befer  to  the  flowcharts  of  Appendix  A. 

A  basic  faailiarity  with  the  HP-87  is  assuaed.  Load  the 
progras  froe  disk  storage,  then  initiate  prograe  operation 
with  the  BOH  key.  The  prograe  will  pronpt  the  user  to  sig¬ 
nify  if  printer  output  is  desired  or  not.  Beply  with  y  or 
B.  The  prograe  code  is  set  for  a  printer  address  of  701. 
Ensure  the  printer  has  this  address,  or  list  the  program 
code  and  change  the  address  in  line  60. 

The  HP-87  allows  a  prograseer  to  set  up  14  predefined 
keys  for  selective  branching.  In  this  prograe,  selections 
are  established  for  each  of  five  coeaonly  used  laeinate  ma¬ 
terials:  T 300/5208,  B (4) /550  5,  AS/3501,  SP1002,  and 
KV49/EPOXI.  Selection  of  any  one  of  these  five  automatical¬ 
ly  inserts  engineering  constants  into  appropriate  variable 
naees.  The  units  of  these  variables  ar?  in  the  SI  system. 
Osers  desiring  to  work  in  English  units  need  only  select  the 


"SI  ->  EHG"  key  prior  to  selecting  a  material.  i  "flag  va¬ 
riable"  is  set  to  unity  for  later  logic  tasting  to  detersine 
whether  or  not  to  apply  a  conversion  factor  to  each  of  the 
input  variables. 

Osers  desiring  to  define  their  own  aaterial  aay  do  sc  by 
selecting  "User  Def'd".  The  program  will  prompt  for  a  ma¬ 
terial  name  and  each  of  the  engineering  constants  before 
continuing. 

The  "Reference"  key  selects  a  display  with  an  abbreviated 
explanation  of  the  program,  a  list  of  the  references,  and 
author  information.  Press  COHriHOB  to  return  to  the  aateri¬ 
al  selection  display. 

Unit  ply  thickness,  ho,  is  then  requested.  Ensure  that 
the  dimension  of  this  parameter  is  the  same  as  that  of  the 
previous  variables. 

Following  this,  and  in  several  subsequent  locations,  the 
user  can  choose  to  display  computed  variables  or  to  bypass 
the  display.  If  a  printer  is  used  with  the  program,  prompts 
for  variables  and  displays  will  not  appear  on  the  CRT,  but 
will  appear  on  the  printer  instead.  Safer  back  and  forth 
between  the  two  to  determine  what  the  next  input  should  be. 
Pith  a  printer  connected,  only  a  question  mark  will  appear 
on  the  CRT.  First  time  users  are  encouraged  to  use  only  the 


CHT  display  (no  printer)  in  order  to  get  used  to  answering 


proapts  and  observing  output.  Run  tine  for  a  complete  ser¬ 
ies  of  calculations  requires  only  a  few  ainutes.  Errors 
should  be  handled  by  stopping  prograa  operation  with  the 
PAUSE  key,  then  reinitialize  with  the  RUH  key.  Effort  has 
been  expended  to  sake  this  prograa  as  user  "friendly**  as 
possible.  Host  proapts  are  self  explanatory. 

Following  optional  displays  of  aodulus  and  coapliance  ma- 
trices  as  well  as  strength  paraaeters,  the  prograa  will  sig¬ 
nify  that  coaputation  of  in-plane  variables  are  now  pro¬ 
gressing.  A  proapt  for  CORB/HO  CORE  laainates  is  then 
displayed.  CORE  laainates  are  those  coaprised  of  a  honey- 
coab  sandwich  section  placed  on  the  centerline  that  contri¬ 
butes  little  to  the  stiffness  of  the  laainate  but  is  used' to 
alter  the  weight  characteristics.  The  thickness  of  the  core 
is  usually  input  as  an  integral  multiple  of  the  ply  thick¬ 
ness.*  Proper  input  of  ply  information  in  sequential  order 
is  very  iaportant  for  flexural  properties  later  ir.  the  pro¬ 
gram.  CORE  and  HO  CORE  materials  both  require  ply  informa¬ 
tion  to  be  entered  from  the  center  ply  outward  in  consecu¬ 
tive  order.  Ho  attempt  has  been  made  to  group  ply 


A  non-integral  auitiple  for  core  thickness,  nc 
tered.  but  the  initial  ply  index  numbs 
have  to  be  entered  as  an  integer.  Use 
pie  rounded  off  to  the  nex*  highe?t 
aay  occur, 
nates  with  a 


may  be  en- 
wiil  still 
nearest  aulti- 

_ _  ... .  _ _ „  _ ,  _  Slight  errors 

but  they  should  be  significant  only  for  lami- 
relativaly  thick  unit  ply  thickness. 


9r,  ti, 
the 
'.umber. 


i 


orientations,  and  the  user  will  be  required  to  enter  each 
ply.  Programming  was  simplified  by  this  expediency,  but  us¬ 
ers  entering  complicated  laminates  or  laminates  with  many 
plies  may  find  it  somewhat  tedious.  Peedback  from  the  prin¬ 
ter  is  immediate,  however,  so  any  mistakes  may  be  discovered 
quickly. 

Following  input  of  CORE/NO  CORE  information,  the  in-plane 
modulus  of  laminates,  Aij,  is  computed  followed  by  its  in¬ 
verse,  aij.  Normalized  versions  of  these  matrices,  Aij*  and 
aij*,  are  computed  as  well  as  the  effective  in-plane  engi¬ 
neering  constants.  A  display  of  these  variables  is  option¬ 
al. 

The  program  then  requests  input  of  loads  (stress  resul¬ 
tants)  N 1 ,  N2,  and  N6,  followed  by  a  request  for  the  orien¬ 
tation  of  the  ply  to  be  examined.  Variables  computed  from 
these  inputs  include  the  in-plane  strain  variables,  on-axis 
strain  and  stress  variables,  max  stress-strain  strength  rat¬ 
ios  (R/R1)#  and  quadratic  strength  ratios  (R/R1).  A  non-op- 
tional  display  lists  results  of  these  calculations.  Follow¬ 
ing  the  display,  the  user  can  enter  a  new  load  and 
orientation  and  repeat  the  calculations  and  display.  If  a 
new  load /orient at ion  is  not  desired,  the  program  proceeds  to 


flexural  computations 


Without  farther  prompts,  the  program  computes  the  flexur¬ 


al  modulus  matrix,  Oij,  followed  by  its  inverse,  dij.  And 
in  a  manner  similar  to  in-plane  computations,  normalized 
parameters,  Dij*  and  dij*,  are  computed  as  well  as  effective 
flexural  moduli,  Eif.  A  display  of  these  parameters  is  op¬ 
tional. 

The  program  now  prompts  for  the  moments.  Ml,  M2,  and  M6. 
An  index  number,  t,  representing  the  sequential  ply  number 
for  the  ply  to  be  examined  (coanting  from  the  ply  closest  to 
the  center  line)  is  requested  along  with  the  orientation  of 
that  ply.  With  these  inputs,  the  program  calculates  curva¬ 
ture,  ki,  on-axis  stress-strain,  max  stress-strain  strength 
ratios  (R  6  R')  ,  and  quadratic  strength  ratios  (R  £  s').  A 
non-optional  display  presents  all  parameters  in  a  format  si¬ 
milar  to  that  for  in-plane  parameters. 

Following  the  display,  the  user  can  enter  a  new  set  of 
moments,  index  t,  and  THETAt  for  a  repear  of  the  calcula¬ 


tions.  Selecting  the  other  option  terminates  program  opera- 


Chapter  III 
SAMPLE  CALCULATIONS 

The  following  pages  are  reproduced  froa  saaple  calcula¬ 
tions  generated  by  the  LAMINATE  prograa.  The  first  prograa 
run  displays  all  optional  displays.  Subsequent  prograa  runs 
delete  aany  of  these  displays  since  T300/5208  was  used  for 
all  calculations  and  aany  of  these  displays  would  not  vary. 
The  ply  designations  for  each  coaposite  are  noted  next  to 
the  title  for  each  run.  The  printer  used  for  these  runs  was 
configured  for  a  laterally  coapressed  node  to  allow  for  re¬ 
port  foraat. 

Many  of  the  proapts  do  not  appear  on  these  print  outs. 
The  Frograa  is  written  to  display  aany  of  the  proapts  only 
on  the  CRT.  The  reader  is  referred  to  the  flowcharts  to  re¬ 
solve  questions  about  prograa  operation  a  specific  points. 


COMPUTATION  OF  DP-PLANE  AND  FLEXUtAL  PROPERTIES  OF  SYNCTRIC 
LAMINATES  ON  THE  HP-87  MICROCOMPUTER 

Selected  cinpesite  net trial  =>  T39I/S2I8 

Inpet  lwinatc  thickness,  he  (SI  tr  END. 


Ex  *  181.  DUE +1  OP 
E»  >  14.3IIE+9I? 
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Es  »  7.17IE+IW 
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On-axis  Qij‘s,  Sij’s,  UiQ’s,  Aij’s,  and  aij's  fir  T3II/52I8 
(Mite  ij»xys). 


(  181.81E+0I9 
•n-axis  Qij  3  ( 

( 


2.8969E+II9  I 

1I.348E+0I9  0 

7.17I9E+II9 


(  5,5249£-4i2 
on-axis  Sij  *  < 

< 


-1.547E-812 

97.I87F-I12 


I 

I 

139.47E-H2 


U1Q  3  76.3A82E+II9 
ICQ  3  8S.732SE+II9 
U3Q  3  19.71I4E+II9 
Wi  3  22.6&/4E+8I9 
USQ  3  26.88ME+889 


(  22.726E+MA 

on-axis  Aij  =  < 

( 

(  44.199E-M9 
on-axis  aij  *  ( 

I 


•3&212E+8I6  I 

1.2933E+0I6  0 

.8962SE+IM 

-12.37E-II9  I 

776.78E-M9  I 

U2I.8E-M9 


Strength  paraneters  in  stress  space,  Fij,  and  strain  space,  Gij. 
for  T3M/S288  (Mote:  14=1*5). 
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The  pngrw  will  mw  empire  in-plane  and  flexiral  pr  inertias  if 
loninatcs  specified  by  the  iser.  Either  COKE  •  r  NO  COKE  1  demotes 
can  be  specified  and  any  nanber  if  plies  defined.  The  »ser  will 
netc  that  this  prifran  takas  a  'brite  fwet'  appriach  by  rty last¬ 
ing  an  entry  fir  each  ply.  This  nethid  alliws  sinpler  prigran 
cide,  ensvres  fewer  errirs,  and  enhances  accirocy. 
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In-plane  cinpitatiins  f illiw. . . 

Hill  the  laninatc  cintain  CORE  ir  NO  CORE? 
C/NC  »>  EM  LINE 


NO  CORE  laninatc:  Enter  V,  the  tital  nmber 
if  plies,  surface  ti  surface. 

n  *  ? 


NO  CORE  laninatc  thickness,  h  =  .  BUS 


Niw  enter  the  ply  iruntatim  if  each  ply: 
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In-plane  par  utters:  Aij's,  aij’s,  engineer inp  c  instants,  and 
nirnalized  Aijt's  and  aijt's  fir  T3II/S2I8  (Mate :  ij*126) 
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Incut  Wads  N1,N2>N6 


Mi  =  i  HZ  *  0  N6  =  0 


Mon,  enter  THETAt ,  the  erientatim  if  the  plf 
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in-plant  strain:  tin  =  20 .835E-M9 
tit  --  -.6282E-M9 
tb  i  =  -.1371E-I19 

Fit  THETAt  »  I  ....  0311/5218) 

tn-axis  strain:  cm  *  2.083SE-H8 
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sto*  *  MI77E+II9 
<»f<qeadratic  *  4  *»» 


Mew  load:  N1  »  l  H2  *  I  *  9 


Mow,  enter  THETAt.  the  orientation  of  the  ply 
to  be  exanined  (in  degrees). 

THETAt  =  ? 


THETAt  *  91 


in-plant  strain:  cli  3  20.835E-I09 
e2i  3  -.6282E-II9 
tit  3  -.1371E-I19 


Fir  THETAt  3  91  ....  (T3II/5288) 

on-axis  strain:  exi  3  -.h2B2E-ll9 
(ft  3  20.831-119 
esi  3  -.i94A£-li8 

m-axis  strtss:  sxt  3  -5.30E+llt 
sfi  3  21.374E+II1 
ssa  3  -UME-II8 


strength  ratiis  0  O' 

nax  strain  1.86E+0I5  i,i5E+816 

mi  stress  1.87E+0I5  i.iSE+016 

qiadratic  1.87E+QI5  1.1X+8I& 

aieraqc  stressts:  sti  3  3.7340E+8IB 
sti*  3  2.2688EHI9 
<=f (qiadratic  R  4  R>)> 

tmnmmummmtmmmmmmmmmmmmmm 


Tht  priqran  will  nw  cmpite  flexiral  properties  fir  tM  laninatt 
defined  darinq  tht  in-plane  definitim.  Pnpertits  calcileted 
inclidc  the  flexiral  stiffness  Mtrix  (#ij),  the  flexiral  cittpli- 
ance  Mtrix  (dij)>  efftctiie  flexiral  laninate  nidili  (Eif),  and 
the  nirMlized  versions  if  stiffness  Mtrices  (tijt  and  dijt). 


Flcxtrel  preperties  Dij,  dij,  tffectm  flexerel  engineering 
censtents,  and  the  nernalizcd  Dijt  and  dij*.  (date:  i j*126) 

(  1.67UE+III  .30176E-II1  4.392JE-I14 

flexarnl  Dij  *  (  3>3il3E~lli  1.I492E*412 

(  7.4MK-H2 

(  5.9957E-DI1  -  5445£-lli  4.1S12E-I13 

flcxeral  dij  *  (  3.8258E+4II  -.4247E“liO 

(  1.3389E+II1 

Elf  »  1.6411E+411 
E2f  *  3.i727E+41l 
*21f  «  .19116 

E6f  a  7.17IIE+II9 

(  1.6I39E+I11  2.8969E+II9  4. 217 IE-413 

flexeral  Dijl  *  (  3.1779E+I10  1.H72E-IIJ 

(  7.17IIE+M9 

(  6.24S6E-912  -.5693E-912  4.3241E-424 

flcxeral  dij*  *  (  3.1S19E-I11  -.4424E-I21 

<  1.3947E-I1I 


Inpet  nenents  H1,H2,H6 


Hi  a  i  H2  *  0  IMt  *  0 


Hew,  enter  the  inde*  neither,  t,  w d  the  erientetiw,  THETAt 
if  the  ply  t«  be  examned. 
t,  THETAt  *  7 


t  »  2  THETAt  «  I 


cirwtwt,  kl  *  5.99S7E-991 
k2  *  -.S46SEHI1 
lb  *  4.1S12E-I13 

Fir  t  *  2  and  THETAt  *  I 

fits' 1  strains,  <lf  *  1.498fC-ll4 
i2f  *-.136ff-IM 
tbf  *  1.9378E-914 

fltx'l  strtssts,  slf  *2.721X^117 
s2f  *  2.9284MIS 
*t  *  7 .  MIX-117 

strcnyth  ratiis  R  I' 

mx  strain  5.SX+991  5.SXH91 

mx  stress  S.S1E+II1  S.S1E+991 

asadratic  4,f5E*991  7.11E+H1 


awrafa  strtssts:  stt  *  l.lMEtNf 
st#»  •  1.79S4E+999 
(■fUaadratic  I  4  I*)) 


Inpat  mm  lMds  1(1,16,114 


Him  MMnts:  HI  *  1  H2  *  9  H4  *  I 


Hm.  tattr  the  indtx  naMar,  t,  and  tht  iritntatim,  THETAt 
if  tht  pis  ti  ht  txanintd. 
t,  THETAt  *  1 


x  *  1  THETAt  *  99 


comtire,  II  «  5.99S7E-I81 
12  »  -.5465E-I81 
16  a  4 . 1S12E-913 


Fer  t  »  i  and  THETAt  *  91 

flei'l  strains,  tif  »  7.4947E-88S 
l2f  *  -.6835-885 
t6f  *  5.1896E-I17 

flti'l  stresses,  sif  *  7.5565+885 
s2f  »  -.il2SE+*07 
s6f  »  -  .5135-115 

strength  rot  its  R  R ' 

mu  strain  5.18E+I81  3.19E+II2 

mu  stress  5.25+881  3.26E+I8R 

qvodratic  5.22E+881  3.I1E+882 

atereqc  stresses:  sti  a  1.2535E+0I9 
ste»  a  7.2139E+489 
(»f (quadratic  R  4  R’>> 


COMPUTATION  OF  I  Iff  LANE  AM  FLEXURAL  PROPERTIES  OF  SYMMETRIC 
LAMINATES  ON  THE  HP-87  MICROCOMPUTER 


Selected  conpositc  Material  »)  T38I/52I8 


[02/902/452 /-452] 


Inpat  laninatc  thickness,  hi  (SI  ir  ENG). 


hi  <  .III12S 


Tht  proqran  Mill  new  cinpite  in-plant  and  fltiiral  priptrtics  if 
laninates  specified  bp  the  iser.  Either  CORE  k  NO  CORE  laninates 
can  be  specified  and  anp  nmber  if  plies  defined.  The  iser  will 
mte  that  this  prtqran  takes  e  'bnte  force*  approach  bp  reqicst- 
inq  an  entrp  fir  each  pip.  This  nethid  allows  sinpler  proqran 
cidc,  enseres  fewer  errors,  and  enhances  acctracp. 


|TTTTTTrrrrr  ttttttttttttttttttttttttt7ttto»oTTTTTTTTTTTTTTTTTTTTTToTTTTTT*«7TTTT1 


In-plane  cmpetatiins  follow. . . 

Hill  the  laninatc  contain  CORE  or  NO  CORE? 
C/NC  =>  EM  LINE 


NO  CORE  laninatc:  Enter  V,  the  total  ninber 
of  plies,  serf ace  ti  surface. 


In-plane  parameters:  Aij's,  aij's,  engineering  c instants,  and 
normalized  Aijl's  and  aijl's  for  T388/S2B8  (Note:  ij=126) 


(  1S.274E+887 
in-plane  Aij  =  ( 

( 

(  71.761E-B1I 
in-plane  aij  =  ( 

( 

Eio  *  6.9676E+818 
E2i  =  6.9676E+I1I 
*21o  *  .29613 
E6o  *  2.6888E+I18 

<  7.6368E+81I 
in-plane  Aijl  =  < 

( 

(  1.43S2E-I11 
in-plane  aijl  =  ( 

( 


4.5215E+I87  1.6869E-88S 
1S.274E+887  4. I288E-8I4 
S.3761E+8I7 

— .2124E-1 08  1.3668E-92I 
71.761E-018  -.53UE-I1V 
1.8601E-888 


2.2647E+81I  8.4343E-883 
7.6368E+818  2.I144E-3I1 
2.688IE+I1I 

-.4249E-811  2.733IE-823 
i .43S2E-I11  -.18621-121 
3.7282E-I11 


Inpot  loads  Ki,N2,N6 


Ml  =  1  M2  *  I  M6  *  9 


Hon,  enter  TUTAt,  the  orientation  of  the  ply 
to  be  examined  (in  degrees) . 

THETAt  =  ? 


THETAt  *  45 


21 


mtnmntmittmtinuimttttumttutunutnttutmistxttmuxttnu 


in-plane  strains  ei#  s7i.76iE*lil 
e2»  *  -.2124E-I08 
(6«  3  1.366E-B2I 

For  THETAt  3  45  ....  (T388/S288) 

in-aiis  strains  exo  *  2.5259E-0I9 
eta  =  2S.25K-816 
tso  *  -.93B8E-888 

on-aiis  stress!  sxe  =  46.655E+881 
Sfi  =  33.45IE+I88 
sso  *  -.6WK+8I2 


strength  ratios 


I' 


mx  strain  1.82E+IB6  1.82E+I86 
nax  stress  1.82E+886  I.I2E+0I6 
qoadratic  6.94E+B85  i . 35E+Q86 


average  stresses;  sto  »  3.47IIE+8I8 
sto*  »  6.7588E+0B8 
(=f (qoadratic  8  &  RM) 

tmmttmmrnmxmmutmmtmmmntummmmtmmmmm 


Hew  load:  Hi  -  i  M2  =  I  N4  *  8 


How,  enter  THETAt,  the  orientation  of  the  ply 
to  be  examined  (in  degrees). 

THETAt  *  ? 


THETAt  *  98 


in-plane  strain:  <!•  3  71.761E-I1I 
<2i  3  2124E-IIB 

tbt  3  1.366E-02Q 

Fir  THETAt  *  91  ....  (T3II/5208) 

on-ms  strain:  exi  3  -.2124E-088 
tfi  3  71.761E-I10 
eso  3  -.iH7E-0i8 

tn-axis  stress:  sxe  3  -3.458E+RI2 
Sfi  3  68.l?iE+flll 
sse  *  -  .7SME-M9 


strength  reties  R  I* 

mx  strain  S.41E+I0S  3.33E+094 

mx  stress  S.87E^0tS  3.ME+086 

qiadratic  5.S2E+I85  2.68E+IW 


average  stresses:  sti  3  2.7612E+0I8 
sti»  3  1.297AE+IB9 
(3f(qeadratic  R  4  R')) 


Men  load:  Ni  3  i  M2  3  I  Hb  =  8 


Moo,  enter  THETAt,  the  irientetiin  if  the  gif 
ti  be  exanined  (in  degrees). 

THETAt  3  ? 


in-plane  5 train  til  *  71.7ME-I11 
<2i  =  -.2124E-II8 
t6t  «  1.344RE-I2D 

Fir  THETAt  *  I  ....  (T3M/52I8) 


in-axis  strain:  txi  *  7.1761EHQ? 

IM  *  >2.i2«-IM 
<$«  *  1.366BE-I2I 

in-«iis  strtst:  ui  =  12.98SE+II2 
syt  >  -i .19IE+9II 
ssi  *  9.7991E-U1 


strength  ratiis  I  I’ 

mi  strain  i .  15E+IQb  l.iSE+IJb 
mx  stress  i.ltE+IM  i.lhEXWl. 

quadratic  1.16W06  1.1XH86 

average  stresses:  sti  =  5.8181£+#I8 
Stl»  S.6542E+0IB 
(=f (qiadratic  R  &  R’)) 


The  pnqran  will  mm  cinpite  flexural  pripertits  fir  the  laninati 
defined  during  the  in-plane  definitions.  Pripertits  calcilattd 
include  the  flcxiral  stiffness  Mtrix  (Dij),  the  flexiral  cmpli- 
ance  Mtrix  (dij),  effective  flexural  laninatc  mdili  (Eif),  and 
the  ntrMlized  versxins  if  stiffness  Mtrices  (Dij*  and  dijt). 


Flexiral  properties  Dij,  dij,  effect  in  flcxirel  engineering 

c instants,  and  the  aerMlized  Bijt  and  dijt. 

(Hite:  ij>126) 

(  7.6842E+II1 

.52164E+M1 

2.6791E+IIB 

flexiral  Dij  3  < 

4.4693E+II1 

2.6791E+III 

( 

8.I651E+II0 

(  1.3241E-I02 

-.13l7E-tl2 

-.3964E-IB2 

flexiral  dij  3  ( 

2.29S9E-9B2 

-.7191-112 

( 

1.2771EH81 

Elf  3  1.1329E+I11 

E2f  3  4.5334E+tlt 

»21f  3  .19177 

E6f  3  1.1747E+B1I 

<  1.1526EH11 

7.B24SE+M9 

4.I187EHI9 

flexiral  Dijt  3  ( 

4.7139E+M 

4.I187E+II9 

( 

1.2I98E+I1I 

(  8.8271E-B12 

-.8718E-I12 

-.2642E-011 

flexiral  di  jt  3  ( 

1.S3I4E-IU 

-.4794E-IU 

( 

8.S132E-I11 

»»»♦•*♦**•♦*»■ 


Inf  it  aments  (11,112,116 


MM  112  3  9  H6  3  I 


Hiu,  enter  the  index  nwber,  t,  and  the  erientation,  TlETAt 
#f  the  olf  ti  be  excnined. 
t,  THETAt  3  ? 


t  =  8  THETAt  *  I 


i  » 


»*>♦•**  t  t  *♦♦♦»**»  t  ♦  *  » 


cereaterc,  kl  »  i.3241E*tt2 
k2  *  13I7E-4K 

k6  »  -.396C-II2 

For  t  «  4  ud  THETAt  *  91 

flex'l  strains,  tlf  *  9.93ISE-N6 
at  >  -.98WE-II6 
ebf  *  -.2973E-II5 

flti'l  stresses,  sif  *  9.99HE+II4 
s2f  »  -.149SE+II6 
sbf  *  -  .2131E+H5 

strength  reties  R  l> 

mx  strain  3.91E+0I2  2.41E+II3 

mx  stress  4.IIE+II2  2.44E+II3 

quadratic  3.89E+9I2  2.HE+9I3 


average  stresses:  ste  *  S.8418E+M8 
sta*  «  3.1496E+II9 
(*f(qoadratic  R  4  R1)) 


tMMHMUMM 


Iipet  new  leads  (91,492,946 

Men  nonents:  Hi  =>  i  (92-0 

Hb  =  9 

Mew,  enter  the  index  nenber,  t, 
it  the  ply  to  be  examned. 
t,  THETAt  =  ? 

«nd  the  orientation,  THETAt, 

t  =  4  THETAt  *  45 

csmtsrt,  kl  ■  1.324iE~H2 
12  «  -.13I7E-II2 
k4  *  -.3M4E-M2 

Fsr  t  *  4  and  THETftt  *  45 

flex'l  strains,  tif  *  6.62l3E~tt6 
«2f  *  -.b5XE-IO() 

«6f  «  -.19KE-I0S 

n«»l  strtSMS,  slF  *  2.6244E+IIS 
*2f  >  1.5815E+II5 
s4f  *  -  .2I53E+II4 

strenqth  rat its  R  V 

nax  strain  9.77E+M2  1.3IMI3 

mx  stress  8.S3E+I82  1.3IE+M3 

qsadratic  6.71E+M2  i .  97E+I43 

mraqe  stresses:  sts  *  1. 1172+119 
StS»  •  2.9S48E+II9 
(sf (qsadratic  R  4  10) 


End  sf  prsqran.  Thank  yst. 


Eiii‘ii'il'1 'lii till! til 'ii'IilJ !•*•!* •!•*• • •liii'i 'i'lil'ii*' ‘ill  * •] 


mt 

In-plane  cinpitations  fellow. . . 

Will  tht  laninate  curtain  CORE  tr  NO  CORE? 
C/NC  *>  E®  LINE 


CORE  laninate:  Enter  V ,  the  tital  lumber  of 
laninate  plies,  end  'sc',  the  total  amber  of  eqii  talent  cere 
plies  that  define  its  thickness,  sirface  te  serfacc. 
n,nc  *  ? 

»  *  i6  nc  s  ib 


OWE  laninate  thickness,  h  *  .814 


Cloutin'  hmN  free  the  center  if  the  loninate,  enter  an 
integer  amber  *ti*  (t  initial)  that  represents  the  censecitiee 
inOei  nmber  ef  the  plf  adjacent  te  the  cere.  (Far  instance, 
if  the  cere  is  e  tettl  ef  the  thick,  cewt  2  epeiMlcnt  plf  lsfcrs 
frrn  the  center  se  that  ti«3). 
ti  «  ? 


ti  *  ? 


Km  enter  the  plf  eriantatien,  in  decrees,  far  each  plf: 
THETA(ti),  THETKti+i)  ...»  THETACn/2) . 
m  ->  CENTO  HOST  PLT  TO  OUTES  SOFAS  PIT  <-  «*  (IHPOTANt) 

THETM  ?  )  «  ? 

-45 

THETA(  II  )  *  1 

-45 


THETA<  il  )  *  ? 

45 

THETA4  12  )  *  ? 

45 

THETAt  13  )  =  ? 

91 

THETA!  14  )  »  ? 

91 

THETA<  IS  )  s  ? 

I 

THETAt  16  )  *  » 


I  ♦  ♦  '  t  **  ♦  M  M  l 


I  ♦  ♦  ♦  0  ♦  l  ♦  0  *  ♦  ♦  4  ♦  * 


I  »  I  *  t  M  *  *  I*  '  t  j 


In-plane  parancters:  Aij's,  aij's,  enqineerinq  constants,  and 
narnalizcd  Aijt’s  and  oijt's  far  T3II/52M  (Nets:  ij»126) 

(  1S.274E*II7 

4.S21SE+M7 

1.6869E-885 

in-plane  Aij  =  ( 

15.274E+8I7 

4.8288E-SB4 

( 

5.3761E+BB7 

(  71.761E-818 

-.2124E-888 

1.3668E-I2I 

in-plane  aij  *  ( 

71.761E-018 

-.S311E-815* 

< 

1.8681E-888 

Ell  >  3.4838E+81S 

E2a  >  3.4838E+S1B 
»2ia  >  .29613 

E6o  >  1.344BE+81B 

(  3.8184E+818 

1.13I4E+I1I 

4.2171E-II3 

in-plane  Aij#  *  ( 

3.8184E+818 

1.II72E-I01 

( 

1.344BE+B18 

<  2.87I4E-I11 

-.8497E-811 

S.4472E-823 

in-plane  aijt  a  ( 

2.87I4E-I11 

-.2124E-421 

( 

7.4484E-I1I. 

Iiipit  loads  Ml, M2,  No 


Ml  *  i  H2  =  •  MS  *  0 


Maut  enter  THETAt ,  the  orientation  of  the  plf 
to  be  etanined  (in  degrees). 

THETAt  =  » 


THETAt  *  45 


in-plane  strain:  cli  3  71 .761E-I1I 
*2*  *  -.21241-IN 
ete  »  1.366C-I2I 

Fir  THETht  *  45  ....  (T3M/S2M) 

•n-exis  strain:  ext  3  2.S2S9E-0I? 

m  >  25.257E-I18 
eso  3  -.93IIE-II8 

•a-axis  stress:  sxt  ■  46.6SSE+M1 
svt  3  33.4SIE+RII 
ssi  3  -  .664BE+M2 


strength  rat its  R  r 

nax  strain  1.I2EHI&  1.02E+II6 

mx  stress  1.02E+IK  1.KE+8I6 

quadratic  &.74E+II5  1.35EMI4 

ateraae  stresses:  stt  3  1.73SIE+II8 
ste>  «  J.37S4E*IW 
(3f(quodrttic  R  4  RM) 


The  srtqran  Mill  mn  cinpitc  flexural  properties  fir  the  laninate 
defined  derma  the  in-plant  definititns.  Pnperties  calctlated 
indtdc  the  flextral  stiffness  Mtrix  the  flextral  ctnpli* 
once  Mtrix  (dij),  effect  its  flextral  laninate  ntdtli  (Eif),  and 
the  nwMlized  versions  if  stiffness  Mtrices  (Dij*  and  diit). 


mtttm 


Fltxiral  priperties  Dij,  dij,  effectiie  flexural  enqmttrinf 
cinstants,  and  the  nirnalized  Dijt  and  dijt.  INite:  ij=126) 


(  4.23A6E+II2 
fltxiral  Dij  *  ( 

( 

(  2.4MSE-8M 
flexural  dij  *  ( 

( 


.7S935E+II2  1.33%£+081 
3.4848E+082  1.3396E+081 
?,5876E+08t 

-  5264E-H3  -.2787E-003 
2.9988E-0S3  - . 34S2E-983 
1 . 9516E-802 


Eif  *  7.6I80E+I1I 
E2f  *  6.2S42E+810 
i21f  *  .21370 

ibf  »  1.783IE+01I 


(  7.93WE+01I 
fltxiral  Dijt  =  { 

( 


1.4238E+I18  2.S117E+90? 
6.5333E+818  2.S1WE+II9 
i.7977Emi 


(  1.3&44E-I11  -.28I8E-I11 
flexural  dijt  =  (  l.SWE-lii 

< 


-.184C-0U 

1841E-I11 

S.6086E~811 


Input  mnents  Mi,H2,H6 


IH--1  »2  *  0  H6  s  j 


Mm,  enter  the  index  nunoer,  t,  and  the  orientatitn,  THETAt, 
if  the  pip  tu  be  exanintd. 
t,  THETAt  =  ? 


t  =  t6  THETAt  *  0 


cereatere,  ki  3  2.4645E-I03 
k2  =  -.52ME-803 
k6  » -.2707E-JI3 


For  t  =  16  and  THETAt  *  I 

flex'l  strains,  eif  3  4.9290E-0I6 
e2f  3  - .  ltS3E— 9IS 
e6f  3  *,5415E-m 


flei'l  stresses,  slf  3  8.931IE+9IS 
s2f  3  3 . 38t9€>4«3 
s6f  >  -.3882E+II4 


strenqtA  reties  I  I’ 


mx  strein 
nax  stress 
qeadratic 


1.6XHI3 

1.68E+II3 

1.62E+II3 


1.68E+I93 

1.68E+II3 

1.83+913 


atcraqc  stresses:  ste  3  6.0642E+OIB 
ste»  3  6.8489E+II8 
(=f (qeadratic  R  6  RM) 


lepet  new  leads  Ai,H2,H6 


Mch  aments:  Hi  *  i  H2  3  I  M6  3  0 


Him,  enter  the  index  neither,  t,  and  the  erientatien,  THETAt, 
if  the  ply  te  he  exanined. 
t,  THETAt  3  7 


t  3  14  THETAt  3  90 


cumtire,  ki  *  2.4645E-Q83 
k2  *  -.S244E-I83 
k6  s  -  .27O7E-O03 

For  t  *  14  and  THETAt  *  91 

flex' 1  strains,  eif  =  4.3129E-406 
*2f  =  -.9214E-OI6 
tbf  *  -.4738E-4I6 

flex'l  stresses,  sif  *  4.1952E+QI4 
s2f  -  - , 1S5IE+I86 
s6f  *  - .  3397E+4I4 

strength  ratios  R  R1 

mx  strain  9.DOE+9IZ  5.S4E+I83 

mx  stress  9.S3E+QI2  S.86E+803 

qiadratic  9.14E+002  4.64E+0I3 

aoeraqe  stresses:  st*  =  3.4284E+8I8 
st *  1.7398E+M9 
(=f (qiadratic  R  4  R')> 

mntttttnttmntmmtmnmmmummmmttmttmumunmnt 


Input  new  loads  H1,M2,H6 


Mew  nonents;  Mi  *  i  M2  *  0  M6  *  8 


New,  enter  the  index  nunber,  t,  and  the  orientation,  THETAt, 
of  the  ply  to  be  exanined. 
t,  THETAt  =  J 


t  =  12  THETAt  =  45 


mumtnttmmitttttmutwttmumntmtmttmtmmtttuttmtnt 

c#r#at#re,  ki  *  2.444SE-I03 
k2  -  -.S2ME-I03 
k6  = -.27I7E-405 

For  t  =  12  and  THETAt  *  45 

flex'1  strains,  elf  -  3.6968E-086 
«2f  *  - .  790IE-096 
ehf  =  -.4B61E-BB4 

flex'l  stresses,  slf  =  1.S861E+8I5 
s2f  *  9.4278E+884 
s6f  =  -.7819E+I05 

strength  raties  R  8' 

nax  strain  2 .11E+I83  2.11E+8B3 

nax  stress  1.93E+8B3  2.11E+883 

quadratic  1.30E+883  2.97E+I83 

aierage  stresses:  sti  =  4.8649E+8I8 
ste'  =  i . 1153E+889 
(=f(q»adratic  8  4  8')) 

mtmmttmnnmwmmmmtttmtnmttimttxtmmtttmututttt 


End  of  prigran.  Thank  ye#. 


Chapter  I? 

COHCLOSI OHS 

The  original  intent  of  this  program  was  to  duplicate  and 
enhance  the  TI-59  prograa  of  reference  1  on  a  popularly 
available  microcomputer.  The  HP-87  is  a  unit  more  than 
likely  to  find  utility  in  many  engineering  organizations  and 
it  was  the  hope  of  the  author  to  contribute  to  the  design  of 
laminates  by  automating  many  tedious  calculations,  present¬ 
ing  them  in  an  easy  to  read  format. 

The  prograa  listing  in  Appendix  B  should  allow  quick  and 
easy  trouble  shooting  of  any  uncorrected  programming  errors. 
Modifications  and  additions  should  be  relatively  simple  to 
program. 

Overall,  this  program  allows  considerably  more  flexibili¬ 
ty  than  the  TI-59  program.  It  should  also  serve  well  as  an 
educational  tool  for  students  and  engineers  desiring  to  re¬ 
view  the  theory  of  symmetric  laminates. 


Appendix  b 
PEOGHAH  CODE 


II  P ASSIZE  16 
21  OPTION  BASE  1 

31  CLEM  I  BISP  8  BISP  ‘D«  yei  Hi sh  ti  how  a  capy  if  the  ivtpit  in  * 

41  BISP  *o  printer?  Y/N* 

SI  HfUT  Al 

61  IF  UPC*  <M)='Y'  THEN  PRINTER  IS  711  ELS  PRINTE1  IS  1 
71  CLEM 

81  PRINT  8  PRINT  ‘CONFUTATION  OF  IN-PLANE  AND  FUXURM.  PROPERTIES  OF  SYMMETRIC 
91  PRINT  *  LAMINATES  ON  THE  HP-87  MICROCOMPUTER 1  I  PRINT  I  PRINT 

110  IISP  I  BISP  I  BISP  8  BISP  "Initialization  if  prvqran  variables.  Standby.' 

III  !  Initialize  variables,  strings,  and  arrays. 

120  INTEGER  I, J,t ,ti,n >nc,FLAC 

130  REAL  EX>EY,MUX>ES,X,](P>Y,YPlS>cli>c2e>e6i,elf>t2f)a6f>slf,s2f)s6f 
140  REAL  M,QXX,QYY,«XY>0SS>SXXtSYY>SXY>SSS»El*>E2o>«21«>EAa,Elf,E2f>*2if/E6^ 

ISO  REAL  AXX,AYY,AXY,ASS,«xi,ayy,axy,ass 

161  REAL  FIX, FYY,FXY)FSS,FX,FY,FXYS,CXX,CYY,CXY, CSS, GX, CY, U1Q.U20, 1130,040,1151 

170  REAL  Mi,N2,N6,Hi,N2,N6,kl,k2,k6,h,he,zs,bs,p,q,a,b,c,nn,iw,SE,Ie,Re,x 

180  REAL  VI, 02, 03,04,010,020,030, 041, THETA, THETAi, BET 

190  REAL  e»,tyetfsi,sxe,sye,ssa,sti, sttpp 

201  REM.  REPSX , REPSXP , REPSY , REPSYP , REPSS .REPSSP 

210  REAL  RSIGX,RSICXP ,RSIGY,RSIGYP ,RSIGS,RSICSP 

221  DIM  A<6,6) ,o(6,6) ,AS(6,6),as(6,6) ,Q(6,6) ,THETA(SI I 

230  DIM  0(6,6) ,d(6,6> ,9S<6,6) 

241  01, 02, V3, 04, 010,020, 030, 041, FLAM* 

250  SE*6894. 97543 
261  FXYS»-.S 
270  FOR  I«i  TO  6 
280  FOR  J«i  TO  i 

290  Q(I,J)*0  8  A(I,J)*0  8  AS(I,J)M  0  a(I,J)«0  (  as(I,J)»0 

301  B(I,J)»0  I  BS(I,J)*I  8  d(I,J)»0  9  ds(I,J)=l 

310  NEXT  J 

321  NEXT  I 

331  GOTO  660 

341  CLEAR  I  FlAG=i 

351  PRINT  8  PRINT  'SI  ->  English  Units  Cinversitn;  A  prtqran  ’flaqJ  has  been* 
361  PRINT  'that  will  key  appreprute  l»9ic  te  einwrt  each  if  the  wives  inpit" 
370  PRINT  ‘aetinatically  in  SI  units  by  prtgran  READ  statements.  If  this* 

380  PRINT  'vptiin  is  selected,  be  svre  tt  inpit  'he',  the  laminate  thickness,' 
390  PRINT  *ie  Enqlish  wits  as  wll.' 

410  PRINT  t  PRINT  'Me#,  press  CONTINUE  t»  retvr*  te  the  CRT  displov  if  the' 

410  PRINT  'material  selectiw.'  9  PRINT  8  PRINT 
420  PAUSE 
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431  GOTO  660 

440  PACESIZE  24  8  CLEM  •  POINT  8  PRINT 

451  PRINT  '  The  LAMINATE  program  is  designed  to  quickly  calculate  in-plane" 
408  PRINT  'and  flexural  strength  and  stiffness  properties  of  selected  * 

470  PRINT  'laminates.  The  program  is  interactive  and  prompts  the  oser  for" 

480  PRINT  'all  required  input  parameters.  Other  pronpts  will  request  the  oser* 
490  PRINT  'to  select  or  skip  a  display  of  various  nudulos  and  compliance  * 

500  PRINT  'matrices  ss  well  as  effectiue  engineering  constants.  Host  other* 

Sll  PRINT  'prompts  are  self  explanatory.* 

520  PRINT 

530  PRINT  'References:  (a)  Tsai,  S.N.,  Hahn,  H.T.,  'Introduction  to  Composite* 

540  PRINT  ‘  Materials',  TECHNOHIC  Publishing  Co.,  d980.' 

550  PRINT  *  (b)  Donaldson,  S.L.,  'Rcuised  Instructions  for  TI-S9* 

560  PRINT  *  Combined  Card/Module  Calculations  for  In-Plane  and* 

570  PRINT  *  Flexural  Properties  of  Symmetric  Laminates',  Air* 

580  PRINT  *  Force  Systems  Command,  June  1982.* 

590  PRINT  8  PRINT  'Program  written  by:' 

600  PRINT  *LCBR  D.R.  Ferrell,  USN* 

611  PRINT  'Department  of  Aeronautics* 

620  PRINT  'Naual  Postgraduate  School" 

630  PRINT  'Monterey,  California  93940', '(September,  1982)' 

640  DISP  ,,, "Press  CONTINUE' 

650  PAUSE 

661  PACESIZE  16  8  CLEM 

670  !  Set  key  labels  for  composite  material  selection  display. 

680  ON  KEY*  1 , 'T300/5208*  GOTO  Cl 
690  ON  RET*  2,'B<4>/5505'  GOTO  C2 
700  ON  RET*  3,'AS/3501*  GOTO  C3 
710  ON  RET*  4,'SP1*02'  GOTO  C4 
720  ON  RET*  S,'RW49/EP0XT*  GOTO  C5 
730  ON  RET*  6, "User  Daf'd*  GOTO  C6 
740  ON  RET*  7, 'Refer cnees*  GOTO  440 
750  ON  RET*  8, 'SI  ->  ENG'  GOTO  340 
760  RET  LABEL 

771  DISP  '  Selection  of  one  of  the  fioe  composites  listed  below  automatically* 
780  DISP  'reads  appropriate  engineering  constants  and  oalues  of  conprcssm,* 
790  DISP  'tensile,  and  shear  strengths  into  variable  names  for  later  conoota-' 
800  DISP  'tions.  All  values  are  in  SI  units.  If  you  anticipate  working  in’ 
910  DISP  'English  units,  press  the  'SI  ->  ENG'  key.  If  you  desire  to  define* 
821  DISP  *a  material  not  listed  here,  select  'User  Def'd'.  The  program  will" 
330  DISP  'then  interactively  request  values  for  Ex,  £y,  vx,  Es,  X,  X',  Y,* 

841  DISP  *Y',  and  S.  Select  'References'  for  an  explanation  of  the  program* 
850  DISP  'and  a  listing  of  references.' 

860  GOTO  860 

870  !  Set  data  pointer  to  read  appropriate  data  string  for  constants. 

880  Cl;  RESTORE  D1  •  GOTO  R» 

390  C2:  RESTORE  D2  R  GOTO  RD 
990  C3:  RESTORE  13  8  GOTO  RD 
910  C4:  RESTORE  D4  8  GOTO  RD 
920  C5:  RESTORE  15  8  GOTO  RD 
930  !  READ  and  DATA  statements  for  the  constants. 

941  RD:  READ  COfrOSITE»,EX,EY,NUX,ES,X,XP,Y,YP,S 


958  IF  FLA04  THEN  GOTO  1168 

961  EX3£X/SE  9  EY=0/S£  8  ES*ES/SE 

971  X»X/SE  I  XP»XP/SE  I  Y«Y/S£  I  YP*YP/SE  8  S=S/SE 

986  GOTO  1161 

996  !  MTA  statenents.  See  pages  19  ond  292,  ref  (a). 

1166  61:  MTA  T300/5208, 181666666666, 16366661166,  .28,7176616616,1516166666,18610 
618 16 ,46661666 ,246661166 ,68666661 

1118  92:  MTA  B(4)/S585 ,264666666666 ,18866666661 , .23 ,8596116661 ,1266166666 ,25116 
86191,61169616,262996196,67066661 

1926  93:  9ATA  AS/3S61, 138696606666, 8966696996, .3, 7108666699, 1447669096, 144796666 
1,81766696,216866906,93889896 

1939  94:  0ATA  SP1092, 38666669968, 8279996966, .26, 4149699666, 1662606669, 616668999, 
31099009,118066966,72699040 

1940  98:  DATA  KV49/EP0XY, 76888086808, 5560888898, .34,2360886680,1406608066,238888 
669,12696996,83069999,34999696 

1988  C6:  CLEAt 

1969  PRINT  *Inpet  none  «f  the  naterial:* 

1876  INPUT  C0HP0SITE4 

1680  PRINT  6  PRINT  ‘Inpat  engineering  emstents: ' 

1679  PRINT  'll,  Ef,  n,  Es‘ 

1166  INPUT  EX,EY,NUX,ES 

1119  PRINT  USING  ‘2(3».3DE),.4D,3D.3DE‘  ;  EX,EY,NUX,ES 

1126  PRINT  I  PRINT  *Inp»t  conpressioe,  tensile,  and  shear  strengths:” 

1134  PRINT  *X,  X',  Y,  Y*  S' 

1148  INPUT  X,XP,Y,YP,S 

1158  PRINT  USING  ,S(3».3BE)“  ;  X,XP,Y,YP,S 

1166  CLEAR  8  PRINT  ‘Selected  cenpesite  naterial  *>  ‘jCOHPOSITEI 

1179  PRINT  P  PRINT  9  PRINT  ‘Inpit  laninete  thictness,  h#  (SI  or  ENG).* 

1186  PRINT  6  PRINT  9  IlfUT  he 
1196  PRINT  t  PRINT  "he  3‘;ho 
1296  MIT  1966 

1218  CLEAR  6  6ISP  6  DISP  *5e  yen  msh  to  see  a  display  if  the  engineering1 
1229  DISP  ‘constants,  and  the  conpressiic,  tensile,  and  shear  contants?' 

1238  DISP  ‘  Y/N‘ 

1249  INPUT  A» 

1289  IF  UPC$  (A9)=*N*  THEN  GOTO  1424 

1269  clear  p  print  6  print  'mntntunwmtmtuttwttmMutmmtm 

1279  PRINT  'Engineering  constants  and  conpressive,  tensile,  and  shear* 

1299  PRINT  ‘constants  for  ‘;C0HP0SITE9 
1298  PRINT 

1364  PRINT  USING  '19X ,5A,30 .3DE‘  ;  'El  *  ‘,EX 
1319  PRINT  USING  ‘18X,SA,X,2D.3DE‘  ;  ’Ey  3  ‘,EI 
1324  PRINT  USING  *i9X,SA,3X,.2D*  j  *«  *  ‘,NUX 
1336  PRINT  USING  *18X,8A,2X,D.3DE‘  ;  ‘Es  3  *,E$ 

1348  PRINT  f  PRINT  USING  ‘11X, 3X*  }  ‘X  *  ‘,X 
1386  PRINT  USING  ‘18X,SA,4D.3DE*  ;  ‘X’  3  \XP 
1368  PRINT  USING  ‘liX,4A,2X,a.3DE‘  ;  *Y  3  *,Y 

1379  PRINT  USING  ‘19X,SA,X,3D.3DE‘  ;  T  3  ‘,Tf 

1380  PRINT  USING  ‘11X,4A,2X,2D.3DE‘  ;  *S  3  *,S 
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1391  print  f  print  'utmmmxtatwtmntmimumtnmuumumttt 
i  print 

1401  JISP  'Use  the  <R0LL>  key  ti  view  CRT  iitpot.’,, ‘Press  CONTINUE’ 

1410  PAUSE 

1420  CLEAR  8  DISP  8  JISP  8  JISP  ’Conpitatiin  of  on-axis  nodilis  and  conpliance' 
1431  DISP  ‘terns.  Standby.' 

1441  N=i/(l-N(JX‘2t<EY/£X» 

1451  BXX*H«EX  8  8YY=NIEY  8  BXY=N*NUXI£Y  8  BSS<1 

1460  SXX*1/EX  8  SYY*i/EY  8  SXY*-<NUX/EX>  8  SSS*1/ES 

1470  !  Cinpote  UiB's,  the  linear  conbinations  if  on-axis  nidili  (Bij,  ij=xvsl 

1480  UiB=<3*QXX+3l8YYt2<QXYt4«SS)/8 

1490  U2U=<QXX-«YY)/2 

1500  U3Q= (BXX+BYY-2tBXY-4<8SS ) /8 

1510  U48=<BXX+8YY+tt«Y-4l8SS)/8 

1520  U5Q=((KX*BYY-2<BXY*4<8SS) /8 

1530  !  Cinpvte  m-axis  A  natrix  fir  single  pi) 

1540  AXX^iXXthi  8  AYY=8YYthi  8  AXY=BXY<bi  8  ASS=BSS*hi 
1550  !  Cinpite  inversion  if  A  natrix  (aij>.  See  pqe  77. 

1560  DET=ASSI(AXXtAYY-AXY*2) 

1570  axx*AYYIASS/JET  8  ayy*AXXIASS/DET  8  axy=-(AXY*ASS/DE7)  8  ass=AXXtAYY/0ET 
1580  OISP  8  DISP  ’Di  yn  wish  ti  see  a  display  if  the  w-axis  Bij's,  Sij's* 

1590  DISP  "lliQ>s,  Aij's,  and  aij's?’ 

1600  DISP  *  Y/N’ 

1610  INPUT  Ai 

1620  IF  UPCI  (AO)=*N*  THEN  SOTO  1880 
1630  PAGESIZE  24  8  CLEAR 

1640  print  8  print  •ttntmnntmtmttmnntmmmutttutmttmttm 
mtmttmmtmv  8  print 

1650  PRINT  ‘On-axis  Bij's,  Sij's,  Uifl's,  Aij's,  and  aij's  for  ‘sCOHPOSITEO 
1660  PRINT  ‘(Note!  ij*xys).’  8  PRINT 

1670  PRflff  8  PRINT  USING  •19X,2A)3D.2DE)2X,D.4DE,6X)A,  ;  ’<  ’,BXX,flXY,’0’ 

1680  PRINT  USING  ‘SX.lBA.liX^B.SDE^A’  ;  ‘m-axis  Bij  =  (’(BYY/O’ 

1690  PRINT  USING  •m.A.aTXjB.ADE’  )  ’<‘,BSS 
1700  PRINT 

1710  PRINT  USING  ‘t9X,2A,D.4DE,2X,2D.3»E,6X)A‘  ;  V  ‘(SXX,SXY,‘0‘ 

1720  PRINT  USING  ’SX,i8A,liX,2D.3DE,6X,A’  ;  ‘on-axis  Sij  =  (*,SYY,’0' 

1730  PRINT  USING  ‘19X,A,27X,3D.2DE’  ;  ’<*,SSS 

1740  PRINT  8  PRINT  USING  ‘13X,6A,2D.4DE’  ;  ’U1Q  *  ’,U1I 

1750  PRINT  USING  ‘13X)6A>2D.4DE‘  ;  ‘U28  *  *,U2B 

1760  PRINT  USING  ‘13X,6A,2D.4DC‘  ;  *U3B  *  ’,U3Q 

1770  PRINT  USING  *13X ,6A,2». 4DE*  ;  *U4B  *  *,U4fl 

1780  PRINT  USING  ‘13X,6A,2B. 4DE*  ;  ’USB  *  *,US8 

1790  PRINT  8  PRINT  USING  *17X,2A,2D,3DE,2X,.5DE,6X,A*  ;  ’<  ‘,AXX,AXY,’ 8’ 

1800  PRINT  USING  “SX^tSA^lX^. 4DE,6X,A*  ;  ’in-axis  Aij  =  C,AYY,’0’ 

1810  PRINT  USING  ‘19X,A,27X,.5l)E’  ;  V.ASS 

1820  PRINT  8  PRINT  USING  ■i9X,2A,2».3DE,2X,3I>.2BE,6X,A*  ;  V  ’,atx,axy,T 
1830  PRINT  USING  •SX,18A,ilX,3D,2DE,6X,A‘  ;  'm-axis  aij  *  (\ayy,‘0* 

1840  PRINT  USING  *l?X ,A?27X,4». BE*  ;  ’<’,ass 

1851  print  8  print  wiwioottwiotittmtwittiiiimiwimiimmttti 
tmmtmtmtmv  t  print 

I860  DISP  ’Use  <RXL>  key  to  vie*  CRT  iotpit.‘,, ‘Press  CONTINUE’ 


48 


1871  PAUSE 

1881  CLEAR  8  PAGESIZE  16  8  SISP  8  DISP  8  MSP 

1891  DISP  *Coap»tation  of  strength  parameters  in  strain  and  stress  space. ' 

190#  DISP  ‘Standby. * 

1911  FXX»1/(X*XP)  8  FYY=i/(Y#YP)  8  FXY=FXYStSfi8  (FXXIFYY)  8  FSS=1/S*2 

192#  FX=i/X-i/XP  8  FY=1/Y-1/YP 

193#  GXX*FXX#QXX*2+2IFXY#8XXI8XY+FYY»8XY*2 

194*  CYY»fXX#8XY»2t2tfXYmXY«JYY+RYIQYY*2 

195*  GXY=fXXI8XX«XYtfXY$(aXX»8YY+8XY*2)+FYYWXYI8YY 

196#  GSS=FSSW)SS*2 

197#  GX=FX*QXX+FYtQXY 

198#  GY=FX«JXY+fYWY 

199#  DISP  8  DISP 

200#  DISP  'Do  yea  vish  to  see  a  display  of  strength  paraneters  in  stress' 

2(11  DISP  'space,  Fij,  and  strain  space,  Gij?' 

2*2#  DISP  '  T/N' 

2(30  INPUT  A# 

2040  IF  UPC*  <A*»=*H-  THEN  GOTO  2251 
2(5#  CLEAR  8  PAGESIZE  24  8  CLEM 

tunmmtmtmv  t  print 

2170  PRINT  ‘Strength  paraneters  in  stress  space,  Fij,  and  strain  space,  Gij.* 
2088  PRINT  'for  'jCONPOSITE*;'  (Note:  ij=*ys).'  8  PRINT  8  PRINT 
2090  PRINT  USING  ,17X,2A,.5DE,2X,2D.3DE,6X,A*  ;  '<  ',FXX,FXY,'#' 

210*  PRINT  USING  'iiX,7A,14X,3D.2DE,6X,A'  ;  'Fij  *  (*,FYY,'#' 

211#  PRINT  USING  *17X,A,27X,3D.2DE'  j  'C,FSS 
212*  PRINT 

2130  PRINT  USING  *12X,5A,2D.3DE'  ;  'Ft  *  *,FX 
2140  PRINT  USING  *12X,SA,2D.3DE'  ;  'Fy  *  ',FY 
2150  PRINT  8  PRINT  8  PRINT 

2161  PRINT  USING  'i7X,2A,2D.3DE,2X,2D.3DE,6X,A'  ;  '(  ',GXX,GXY,'#'  ' 

2170  PRINT  USING  'iiX,7A,i4X,2D.3DE,6X,A'  j  'Gij  =  (•,GYY,'I' 

2180  PRINT  USING  '17X,A,27X,2D.3DE'  j  '<',GSS 
2190  PRINT 

2200  PRINT  USING  ,12X,5A,D.4DE'  ;  *Gx  *  ',GX 
2210  PRINT  USING  *12X,SA,D.4DE'  ;  'Gy  =  \GY 

2220  print  8  print  •tmtttttmmmmttttutmtmttmmmttumuttm 

2230  DISP  'Use  the  (ROLL)  key  to  view  CRT  igtp«t.',,'P.tS5  CONTINUE* 

2240  PAUSE 

2250  PAGESIZE  16  0  CLEM  8  PRINT 

2260  print  8  print  •ttmmumnmttmtwnmtmnutttnnmtmmttt 
tuutnmttmmv 

227*  print  ‘titmtmitmtmimmtimmtmmmmmmmtmimmi 
mttmmv  0  prut 

2280  PRINT  ‘The  progran  will  now  conpote  in-plane  and  flexural  properties  of* 
2290  PRINT  'laninates  specified  by  the  user,  Either  CORE  or  NO  CORE  laninates’ 
230#  PRINT  'can  be  specified  and  any  nvnber  of  plies  defined.  The  user  will* 
231#  PRINT  'note  that  this  progran  takes  a  'brote  force'  approach  by  reqoest-” 
232*  PRINT  *ing  an  entry  for  each  ply.  This  nethed  allows  simpler  progran  * 

233#  PRINT  'code,  ensores  fewer  errors,  and  enhances  accuracy.’ 
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2341  PtINT  8  PRINT  *1 


^ZZTm'mmuutn  *  m  ***** 

2361  PRINT  ‘In-plane  ciaputatiins  filliu...' 

2371  PRINT  I  PRINT  'Will  the  laninate  cm  tain  CORE  ir  NO  CORE?' 

2381  PRINT  1  C/NC  «>  END  LINE*  9  PRINT  f  PRINT 
2391  INPUT  M 

2401  IF  UPCt  (At)-*'  THEN  GOTO  2S60 

2411  CLEAR  8  PRINT  f  PRINT  'NO  CORE  laninate:  Enter  V,  the  total  nmber' 

2421  PRINT  'if  plies,  airfare  t»  sirface.' 

2431  PRINT  '  n  *  ?' 

2441  INPUT  n 

24SI  PRINT  8  PRINT  *n  =';n  I  PRINT  8  PRINT 
2460  h=nihi  !  Thickness  if  NO  CORE  laninate . 

2470  PRINT  8  PRINT  8  PRINT  'NO  CORE  laninate  thickness,  h  *';h 
2481  PRINT  8  PRINT  8  PRINT  *Nw  enter  the  pip  irientatim  if  each  ply s“ 

2490  PRINT  'THETA(l),  THETA42),  ,  TKTA(n/2)' 

2500  PRINT  'tit  ->  FROM  CENTER  TO  OUTER  SURFACE  (-  ttt  (1HP0RTANT!)" 

2510  PRINT  8  PRINT  8  PRINT 

2520  ti=i 

2S30  GOSUB  THETA 

2540  zs»t  !  Vilune  fractiin  if  core. 

2550  GOTO  2800 

2560  CLEAR  8  PRINT  8  PRINT  'CORE  lwinate:  Enter  'n\  the  tital  nmber  of* 

257ti  PRINT  'laninate  plies,  and  'nc*,  the  tital  nmber  if  egoiwlent  core* 

2580  PRINT  'plies  that  define  its  thickness,  sirface  ti  ssrfece.* 

2590  PRINT  '  n,nc  =  »• 

2600  INPUT  n,nc 

2610  PRINT  8  PRINT  'n  *';n-,'  nc  =';nc 

2620  WAIT  1000  8  h=(n+nc)thi  !  Thickness  if  CORE  laninate. 

2630  PRINT  8  PRINT  8  PRINT  ‘CORE  laninate  thickness,  h  >*;h 
2640  WAIT  1000  8  CLEM 
2650  PRINT  8  PRINT 

2660  PRINT  *Comtinq  utuard  frtn  the  center  if  the  laninate,  enter  an  ' 

2670  PRINT  'integer  amber  'ti*  <t  initial)  that  represents  the  consecitm" 
2680  PRINT  'index  nmber  of  the  ply  adjacent  ti  the  core,  (Fir  instance,' 

2690  PRINT  'if  the  core  is  a  tital  if  4hi  thick,  cunt  2  eqiiialent  ply  layers' 
2700  PRINT  ‘frin  the  center  si  that  u*3).‘ 

2710  PRINT  *  ti  «  ?• 

2720  INPUT  ti 

2730  PRINT  8  PRINT  8  PRINT  'ti  *';ti 
2740  WAIT  1000  8  PRINT 

2750  PRINT  8  PRINT  'Now  enter  the  ply  irientation,  in  degrees,  fir  each  ply:* 
2760  PRINT  'THETA(ti),  THETA(ti+i)  ...  ,  THETA(n/2).* 

2770  PRINT  'ttt  -->  CENTER  HOST  PLY  TO  OUTER  SURFACE  PLY  <-  ttt  (IMPORTANT)' 
2780  GOSUB  THETA 

2790  :s*ncthi/h  !  Vilone  fraction  if  core. 

2800  CLEAR  8  MSP  8  OISP  8  MSP  'Cmpitatun  of  in-plane  properties.  Stinobi.' 
2810  !  Cmpotation  if  in-plane  nidulis  if  laninates. 

2820  A<l,l)*UiQt(nthi)mtU2W2tU3Q 
2830  A(2,2)=UiQI(ntho)-VttU2(H2*J38 
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2841  A(l,2)*U4QS(nlhe)-V2tU3fl 
2850  A(6,6)=U5«(n*he>-V2*U3Q 

2861  A<i,6)=V3*U2G/2+V4«U3Q 

2870  A(2,6)sV3IU2G/2-V4tU30 

2881  !  Conpvtation  of  in-plane  conplianct  of  laninates  (inversion  of  A( i , j) ) 

2890  DETsA(i,l)tA(2,2)tA(6>6)+2tA(l,2)tA(2,6)tA(i,6)-A(2,2)tA(i,6)A2-A(6,6)SA(l| 
2)A2-A(1,1)«A<2,6)A2 

2900  a(i,i)*<A<2,2)tA(6,6)-A(2,6>A2>/8ET 

2910  a(2,2)*(A<l,t>tA<6,6>-A(i,6)A2>/0ET 

2920  a(i,2)*(A(l,6)»A<2,6>-A<i,2)IA<6,6))/0ET 

2930  a(6,6)*(A<l,l)IA<2,2>-A<l,2>A2)/0ET 

2940  a(l ,6)*(A(1,2) «A(2,6)-A(2,2)«A( 1 ,6) ) /DET 

2950  a(2,6)*(A(l,2)tA(l,6)-AU,l)tA(2,6))/DET 

2960  !  Conpotaiion  of  ntrnalited  in-plane  nedvlus  and  conpliance 

2970  FOR  1*1  TO  6 

2980  FOR  J*1  TO  6 

2990  AS<I,J)*AU,JVI 

3000  as(I,J)*a(I,J)lh 

3010  NEXT  1 

3020  NEXT  I 

3030  !  Coapvtation  of  effective  in-plane  engineering  constants. 

3040  Eie*l/(a<l,i)*h)  !  Longitvdinal  nodslis 

3050  E2a*l/(a<2,2)tli)  !  Transverse  nod  ulus 

3061  E6o*l/(a(6,6)Hi)  !  Shear  nedvlis 

3070  v2io*-(a(l,2)/a(l,l))  !  Poisson's  ratio 

3080  819  t  OISP  I  OISP  ‘On  poo  wish  to  see  a  display  of  effective  in-plane  " 
3090  OISP  ‘engineer inq  constants,  the  in-plane  stiffness  natrix  (Aij)* 

3110  OISP  ‘and  its  inverse  (aij),  and  the  nornalized  versions  of  these* 

3110  OISP  ‘Matrices  (Ai jl  and  aijt)?‘ 

3120  OISP  ‘  Y/N‘ 

3130  INPUT  AS 

3140  IF  UPCS  (AS)**N‘  THEN  GOTO  3430 
3150  PAGESIZE  24  8  CLEAR  8  PRINT 

3160  print  (  print  •tttmtmmttmmtutmtttttttmtttwtmtmtmtm 
muvnutttmtttv  i  print 

3170  PRINT  ‘In-plane  paroneters!  Aij's,  aij's,  engineering  constants,  and’ 

3180  PRINT  ‘nornalized  Aijl's  and  aijt's  for  *;C0HP0SITE0}*  (Note:  ij*126)* 

3190  PRINT  8  PRINT 

3200  PRINT  USING  ‘19X,2A,2O.3OE,2X,l.4OE,2X,0.4DE‘  j  ‘1  ‘,A(l,i),A(l,2),A(l,6) 
3211  PRINT  USING  ■4X,i7A,i3X,2».31J€,2XJD.40€‘  ;  ‘in-plane  Aij  =  <‘,A<2,2),A<2,6I 

3220  PRINT  USING  ‘19X,A,27X,0.4DE‘  j  V,A<6,6> 

3230  PRINT 

3240  PRINT  USING  ‘19X,2A,2D.31>£,2X,D.«DE,2X,I>.4DE‘  j  ‘(  ‘,o(l,l),a(i,2),a(l,6) 
3250  PRINT  USING  ‘4X,17A,i3X,20.30E,2X,0.40E‘  j  ‘in-plane  aij  =  C,a(2,2),a(2,6l 

3260  PRINT  USING  '19X,A,27X,0.«E’  ;  ‘(‘,a(6,6) 

3270  PRINT 

3280  PRINT  USING  *13X,6A,D.4K*  ;  ‘Ein  *  *,Eio 

3290  PRINT  USING  ’i3X,6A,0.48E‘  ;  *E2o  *  *,E2o 

3300  PRINT  USING  ‘12X,7A,.50‘  ;  *v2to  *  ‘,»21o 
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3318  HINT  USING  •i3X.6M.4DE'  J  *E6e  3  *,E6» 

3321  PRINT 

3331  PRINT  USING  ,l?X,2A,D.4D€f2X,D, 4De>2X,D.4DE*  ;  *(  \AS(i,i),AS(i,2>,AS(i,6> 

3341  PRINT  USING  "3X , i7A,  14X ,B. 4DE ,2X >* . 4DE"  ;  *in-plone  Aij*  *  ( • , AS< 2,2) .AS<2, 
6) 

33SI  PRINT  USING  *19X,A,27X,D.4DE*  ;  *(*,AS<6,6) 

3361  PRIM 

3370  PRINT  USING  '19X,2A,D.4DE,2X,D.4DE,2X,D.4DE*  ;  *<  * ,as(i,i> ,os(i ,2) ,as(  1 

3381  PRINT  USING  *3X,i7A,l4X,D.4DE,2X,D.4IE*  ;  *in-plane  aij*  •  <\os(2,2),os<2, 
6) 

3398  PRINT  USING  •19X,A,27X>».4I>E-  j  *<*,as<6,6) 

34iJ*SS*^S*<Su!>  ktl*!!7*!**  CRT  o  vtp  it.  \, ‘Press  CONTINUE* 

3428  PAUSE 

3438  CLEAR  8  PACESI2E  16 

3448  DISP  'It  util  aw  be  necessary  ti  inpit  the  loads,  Ni,  N2,  aid  N6.* 

3458  DISP  'Outsider  ttw  eptions:* 

3468  DISP  *  Disc  (a):  Inpat  selective  snit  leads  to  dettrnine* 

3478  DISP  *  nai mm  leading  alleuahle.* 

3488  DISP  '  Casa  (b):  Inpet  icteal  loading  to  detcrnine  strain* 

3498  DISP  '  invariants  and  strength  ratios.* 

3508  DISP  '(The  pregrut  operates  the  sane  in  either  case.  It  is  up  to  the  * 

3518  DISP  'osar  to  interpret  the  evtpvt  properly).' 

3528  PRINT  I  PRINT  8  PRINT  'Inpvt  leads  N1,N2,N6*  8  PRINT  8  PRINT 
3538  INPUT  N1,N2,N4 

3540  PRINT  8  PRINT  *N1  **;Ni;*  N2  **;N2;*  N6  **;N6 
3558  MIT  1888 

3568  !  Conpote  the  in-plane  strain  variables,  eie,  e2o,  t6o.  (T4.2,  p . 121 ) 

3578  ele*a<i,i)*Nl*a<l,2)tN2*a(i,6)IN4 

3588  e2o«a(i,2)*Nl*a(2,2)*N2+a<2,6)*N6 

3598  c6o*a(l,6)tNi+a<2,6)fN2+a<6,6)SN6 

3608  PRINT  8  PRINT  8  PRINT  8  CLEAR 

3618  PRINT  8  PRINT  *Neu,  enter  THETAt,  the  orientation  of  the  ply* 

3628  PRINT  *to  be  exanined  (in  degrees).' 

3638  PRINT  *  THETAt  =  ?*  8  PRINT  I  PRINT 
ii4l  Ttwirr  THETAd 

3658  PRINT  8  PRINT  'THETAt  =*;THETAd 

3668  dAIT  1888 

3678  THETAt*T8CTAd8PI  /188 

3688  CLEAR  8  DISP  8  DISP  8  DISP  'Conpvtation  of  in-plane  properties.  Standby* 
3698  !  Conpote  the  on-axis  strain  variable,  ext,  eye,  eso,  for  the  given  ply. 
3780  !  T2.6,  p.53 

3718  p*(eie*e2o)/2  I  q*(elo-e2o)/2  8  r*e6e/2 

3728  m«C8S  (2ITHETAt)  8  nn*SIN  (2ITHETAt) 

3738  exo*p+q*nn+r*nn 

3748  eyo«p-glM-rlnn 

3758  eso*-(28glnn)+2lrtnn 

3768  !  Conpote  the  tn-axis  stress  variables,  sxo,  syo,  sso,  for  the  given  ply. 


r* 
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3771  «  Ti.6,  p.14 
3781  sxe*ttXleieo«XYIeM 

3791  spe«GXYIexe4ffrY*«pd 

3811  sse«8SSIese 

3811  G08U8  SSSt  !  Co  t«hr  t«  conpotc  the  stress-strain  strenqth  rat  us,  R4R' 
3828  C0SU8  OUAI  !  Co  sibr  ti  conpote  the  qoadrotic  strenqth  ratiis. 

3838  stoMq/h  •  stopp»*qp/h 
3841  CUM  f  PACES12E  24  »  CUM  I  PI  1*1 
385IP«MT8P8m  #  ^ 

3868  PRINT  USDS  •18X,23A,2».3D€«  ;  'in-plane  strain:  ela  *  Veil 
3878  PUNT  USING  ,27X1iAtD . 4P£*  j  *e2e  *  \e2o 
3888  PRINT  USING  •RTX.M.D.ODE'  }  'e6e  *  \e6o 

3898  PRINT  8  PRINT  for  THETAt  »';THETAd}*  ....  (*;COHPOSITEi;*)'  I  PRIil 
3988  PRINT  USING  -28X,22A,D.40E’  }  'on-axis  strain:  tie  *  ',exe 
3918  PRINT  USING  *36X,6A,2i)  3JE*  }  'ope  *  %epe 
3928  PRINT  USING  '36X^,0.408'  ;  'eso  =  Vs* 

3931  PRINT 


'on-axis  stress:  sxo  =  \sxe 


.5»» 

■,  'sse  *  \sse 

•strenqth  ratios',  I*,  ■R>l 


3940  PRINT  USING  '20X,22A,2D.3DE* 

3958  PRINT  USING  '36X,6A,2D.3RE'  ;  'spa  * 

3968  PRINT  USING  *36X,6A,D.4K' 

3978  PRINT  8  PRINT 

3988  PRINT  USING  *ltX f ISA r8X,A,l 1X,2A*  , 

3998  PRINT 

4888  PRINT  USING  •15X,10A,4X,D.D1)E,3X,8.8D£*  ;  *nax  strain  *,Re,Rep 
4818  PRINT  USING  '1SX,10A,4X,D.DDE,3X,8.»OE'  ;  *nai  stress  Ms,Rsp 
4828  PRINT  USING  *lSX,iOA,4X,D.DDE,3X,D.D8E'  j  'qeadratic  *>Rq>*qp 
4838  PRINT 

4848  PRINT  USING  “2»X,24lh,B. 4DE*  ;  'aoeraqe  stresses:  sta  3  ',sto 
4858  PRINT  USING  '37X,7A,D.4DE'  }  'iti*  3  *,stopp 
4868  PRINT  USING  181, 22A*  }  '<=f(q»adratic  R  4  R*))' 

4878  print  i  print  ■mtttmwwtiiwmttttmttttmmmmmwtmttt 
tttmtmnttmuv «  print 

4888  MSP  *llse  (ROLL)  kep  te  eieu  CRT  oetpet,', /Press  CONTINUE* 

4890  PAUSE 

4180  CLEAR  8  PAGESIZE  16 

4118  MSP  *t>e  pan  wish  to  conpote  in-plane  properties  for* 

4128  DISP  *a  new  load  and  THETAt V 
4138  MSP  *  Y/N' 

4148  INPUT  M 

41S8  IF  UPC*  (AI)**N*  THEN  GOTO  4281 
4168  MSP  8  DISP  'Enter  new  Ni,tt,N6.' 

4178  INPUT  Nl,N2,N4 

4188  PRINT  8  PRINT  I  PRINT  'New  load:  Ni  *';N1/  H2  *-;N2;*  N6  =*jN6 
4198  UNIT  1888  8  GOTO  3S68 
4288  CLEAR  8  PRINT 

4218  print  8  print  •ro*mt»mmummumm*inroi«mntmmim** 

tmmmmmuu  * 

4228  print  •tnntnnnxutnuuwttnmninuutnnutnnututttunttt 
— - 1 PRINT 


4238  PRINT  ’The  preqro*  will  now  conpote  flexural  properties  for  the  laninatt* 


t 

7# 
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4241  PUNT  'defined  firing  the  in-plene  definition.  Properties  coicoiated* 
42S0  PRINT  'include  the  flcioral  stiffness  Mtrix  (Dij),  the  flexortl  conpli-' 
4261  PRINT  'once  Mtrix  (dij),  effect  in  flexors!  lsninote  nodoli  (Eif),  and* 
4271  PRINT  'the  normlixed  versions  of  stiffness  Mtriccs  (Si jt  and  dijl).' 


4311  PAUSE 

4321  CLEAR  »  DISP  9  DISP  I  DISP 

4331  MSP  'Conpotetion  of  flexorol  properties.  Standby." 

4341  hs»(l-is*3)lh43/12 

43SI  !  Conpote  flexortl  nodolos  (TS.3,  p.t80> 

4361  D(  1  ,1  )»WQ*s*V19«U29*V2D«U38 
4371  9<2,2)*lll|lhsH/19tU2W»tU3Q 

4380  9<  1 , 2)*U4Qths-V29IU3Q 

4391  9(6,6>=U5Q»is-V28IU38 

4410  D(l,6)>U2QtV3M24U4MU38 

4411  9(2,6)»«2Q«V»/2-M48*U3e 

4420  !  Conpote  flexorol  coup lionet  < interne  of  Si j) 

4430  DET3D(1, Ot9(2, 2)19(6, 6)+2t9(l, 2)09(2, 6)09(1, 6)~9(2, 2)09(1, 6)A2-D(6,6)tD 

<1,2)A2-0(1 , 1)09(2, 6>A2 

4440  d(l,l)3<9(2, 2)09(6, 6)-0(2,6)A2)/9ET 

4450  d(2, 2)  *(9(1, 1)09(6, 6HH1,6)A2)/9ET 

4460  d(i,2)»(9(l,  6)09(2, 6)-9(l,  2)00(6, 6))/9ET 

4470  d(6,6)=(9(i, 1)19(2, 2)-9<l,2)A2)/9ET 

4480  d(l,6)-(9(l , 2)00(2, 6)-9(2, 2)09(1 ,6) )/D£T 

4490  d(2,6)«(9(i, 2)09(1, 6)-9(l, 1)09(2, 6))/0ET 

4500  !  Conpote  norMlized  versions  of  91  j  and  dij. 

4510  FOR  M  TO  6 
4520  FOR  J»1  TO  i 
4530  98(1 ,J)30(I , J)012/h *3 

4540  ds(I,J)*d(I,J)»43/12 

4559  NEXT  J 

4560  NEXT  I 

4570  !  Conpote  effective  flexorel  lsninote  nodoli. 

4580  Elf»12/<hA30d(l,l)> 

4590  E2f=12/(hA30d(2,2)) 

4600  E6f=12/(hA30i(6,6)> 

4610  v21f*-(d(i,2)/d(i,i)) 

4620  DISP  0  DISP 

4630  DISP  *9o  yov  wish  to  see  a  display  of  flexortl  nodolos  and  conplionct* 

4640  DISP  'Mtriccs  (9ij  and  dij),  effective  flexorol  Imitate  nodoli,  and* 

4650  DISP  'the  norMlized  versions  of  the  nodolos  and  conplisnce  Mtriccs)* 

4660  DISP  '  Y/N* 

4670  INPUT  At 

4680  IF  UPC*  <At)«*N'  THEN  SOTO  4970 
4690  PA6ESIZE  24  I  CLEM  I  PRINT 

4700  print  9  print  t$smnt$mmuutmmttmmntttmttm»wmuu 
mtmxttmmrat' »  print 
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4711  PUNT  'Fleioral  properties  Dij,  dij,  effect i»e  flexoral  enqineerinq* 

4721  PRINT  ’constant*,  and  the  normalized  Dij*  and  dij*.  (Note:  ij»i2&)’ 

4731  PRINT 

4741  PRINT  USING  a19X>2A>».4De,2X,.S0C>2X>D.4DEa  ;  '<  \D(1,1>,D<1,2>,D<1,6> 

4751  PRINT  USING  *4X,16A,14X>D.40C,2X>1>.4K>  )  'flexors!  Dij  *  (', 0(2,2), 9(2, 6) 
4761  PRINT  USING  '19X,A,27X,D.48E'  ;  a(a,D(6,6) 

4771  PRINT 

4781  PRINT  USING  'i9X,2A,D.4DE,2X,D.4DE,2X,D.4DE'  ;  ’(  ',d(i,i),d<i,2),d(l,6> 
4791  PRINT  USING  '4X,16A,14X,D,4RE,2X,D.4DE*  ;  ’flexoral  dij  *  C,d(2,2),d(2,6l 
4891  PRINT  USING  't9X,A,27X,D.49E*  ;  '(',d(6,6) 

4811  PRINT 

4821  PRINT  USING  *13X,6A,D.4DE'  ;  "Elf  *  ',Elf 
4831  PRINT  USING  '13X,7A,D.4RE'  }  aE2f  *  a,E2f 
4849  PRINT  USING  *12X,9A,,59'  ;  a»21f  *  a,*21f 
4859  PRINT  USING  '13X,6A,D.41E'  ;  *E6f  -  a,E6f 
4869  PRINT 

4879  PRINT  USING  *i9X,2A,D.4DE,2X,9.40E,2X,D.4DE'  i  *<  a>DS(i,l),DS(l)2)>DS(l,6) 

4889  PRINT  USING  '3X,17A,14X,D.49E,2X,8.4DE'  j  ’fletirel  Dijt  *  C,9S<2,2>,DS(2, 
6) 

4899  PRINT  USING  '19X,A,27X,D.49E'  ;  'C,DS(6,6J 
4999  PRINT 

4911  PRINT  USING  '19X,2A,D.4DE,2X,D.4D£,2X,D.4DE'  }  a<  *,ds(i,i),ds(i,2),ds(l,6> 

4929  PRINT  USING  a3X,i?A,i4XJD.4IE,2XJD.4Ka  }  afleior«l  dijt  *  r,<M2,2Ms<2, 
6) 

4939  PRINT  USING J19X,A,27X,P.4BE*  ^  ^  ^  ^  ^ 

4959  DISP  'Use  (ROLL)  key  to  view  CRT  oitpot,*,, "Press  CONTINUE* 

4969  PAUSE 

4979  PAGES1ZE  16  8  CLEAR 

4989  DISP  a!t  uill  now  be  necessary  to  inpot  moments  HI,  H2,  and  N6.* 

4991  DISP  ’Consider  too  options.' 

5999  DISP  *  Case  (a):  Inpot  selective  onit  moments  to  determine  the* 

5119  DISP  *  max  loadinq  allowable.* 

5929  DISP  '  Case  (b):  Inpot  octoel  moments  to  determine  strenqtn  ratios* 

5939  DISP  *  and  c«r satires.’ 

5949  DISP  '(The  proqran  operates  the  same  in  either  case.  It  is  vp  to  the* 

5959  DISP  *oser  to  properly  interpret  the  reselts).* 

5969  PRINT  8  PRINT  'Inpot  moments  H1,H2,H6'  I  PRINT  8  PRINT 
5979  INPUT  H1,M2,H4 

5989  PRINT  I  PRINT  'HI  *’;Hi;'  H2  *';H2;*  H6  »';H6 

5999  WAIT  1999 

5199  CLEM  8  PRINT  8  PRINT 

5119  PRINT  'Now,  enter  the  index  nomber,  t,  and  the  orientation,  THETAt,' 

S129  PRINT  'of  the  ply  to  be  examined. ' 

5139  PRINT'  t,  THETAt  »?' I  PRINT  I  PRINT  8  PRINT 
5149  INPUT  t,THET4d 

5159  PRINT  8  PRINT 't  »';tj*  THETAt  *';T!€TAd 
5169  UAn  1999 
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5171  TNETAt*TH£TAd»I  /18I 

5181  CLEM  I  DIM  I  DIM  I  DIM  *  Standby* 

5191  !  Cmpitatiia  if  cirntirts,  ki.  (T5.2,  p .  1741 
5201  ki«dd,i)»Nl*d(i,2)tie+dd,6)*H6 
5211  k2»dd,2)0lltH(2,2)IH2*d(2,6)0M> 

5221  k6>dd,6)0Hl+d(2,6)tft2+d<6,6)IM6 

5231  tli*klttthi  !  actially  *lf;  celled  tit  fir  sibr  QUA1 

5241  t2t*k2itOM 

5251  c6i>k6tt4bi 

5261  elf*tli 

5271  t2f*t2i 

5281  c6fst6t 

S29I  !  CMMtt  iff-axis  ciffltMct  fir  stress  calcilatim. 

5311  VWuOCOS  <2ITHETAt> 

5311  V^itCOS  (4tTlCTAt) 

5320  V3*ilSM  <2*TRETAtl 

5330  V4*hilSIN  (4STHETM) 

5340  A(l,l)*U18ftHVl«U2Q«V2SU3Q 

5350  Qd,i)*A(i,i)Al 

5360  A(2,2>»UiQ0bH>10U28»V2*U3Q 

5370  Q(2,2)-A(2,2)/hi 

5380 

5390  Q(l,2)*A(l,2)/hi 

5410  A(6,6)>U5(Xhi-V2«U3l 

5411  Q<6,6)*A(6,6)/hl 

S420  Ad,6>=U2QW3/2»V4IU3l 

5430  Qd,6)«AU,6)/!ii 

5440  A(2,6)«U2CtV3/2-V4*U30 

5450  Q<2,6)«A<2,6)/hi 

S460  slfrQ(t,l)tetf+Qd)2)te2f+Q(l,b)*e6l 

5470  s2f»Q( 1 ,2)0(14 +8(2 ,2)tc2f +Q(2 ,6)ta6f 

5480  s6fa8d,6)»t6f*8<2,6)ft2f*Q<6,4>*a6f 

5490  !  Cmpute  on-axis  strain  fir  R/R*  calcilatim. 

5500  !  (T2.6,  p.53) 

S510  p»<eli*2e)/2  i  q*<elo-t2i)/2  0  r*«6i/2 

5520  mpCOS  (2lTHETAt)  0  nn«SIN  (2tTHETAt) 

5530  txi^qlMHrkmi 

5540  eyi^-qlnrrOnn 

5550  esi»-(2lq*nn)+2<r*nn 

5560  !  Ciflfitt  m-aiis  stress  fir  R/t*  calcilatim, 

5570  sxHKXfdo+OXYttyt 
5580  3>i=QXY0«ii+0YY«aM 
5590  ssi*gSSttsi 
5600  GOSUB  SSSR 
5610  COSUB  QUAD 

5620  sti*6IRq/h*2  I  stipp*6»Rqp/b‘2 
5630  CLEM  0  PA6ESIZE  24  0  CLEM 

5650  PRINT  USING  *11X)16A)D.4DE*  ;  'cimtife,  kl  *  *,kl 
5660  PRINT  USINC  *21X,SA,D. 4D€*  ;  *k2  «  \k2 
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5471  PRINT  IISINC  a21X,5A,B.4DEa  j  ak4  *  a,k4 

5481  PRINT  «  PRINT  fir  t  *a;t;aond  TNETAY  *a;THETAd  I  PRINT 

5491  PRINT  USING  a4X,22A,D.4BEa  ;  Vlei'l  strains,  elf  >  \eit 

5711  PRINT  USING  a28X,4A,B.4BEa  ;  at2f  >  \e2f 

S7it  PRINT  USING  a28X,4A,B.4Ka  ;  *e4f  *  a,«4f 

5721  PRINT 

5738  PRINT  USING  'SX^.IMDE'  i  'fltx'l  stresses,  slf  =  a,slf 
5741  PRINT  USING  *2IX,4A,0.49ea  ;  as2f  *  \s2f 
5751  PRINT  USING  a28X,4A,B.4BEa  ;  as4f  *  a,s4f 
5741  PRINT 

S77f  PRINT  USING  alSA,8X,A,ilX,2Aa  ;  ‘stnnqth  r*tiisa,aRa,aR,a 
5781  PRINT 

5791  PRINT  USING  a5X,18A,4X,B.BBE,3X,8.BD£a  ;  aMx  strain  \Re,Rep 
5811  PRINT  USING  aSX,18A,41,D.BBE,3X,B.D8Ea  ;  aMX  stress  a,ts,Rsp 
5811  PRINT  USING  a5X,18A,4X,B.DBE,3X,D.BB£a  ;  'qiodratic  a,Rq,Rqp 
5821  PRINT 

5831  PRINT  USING  a28X,24A,B.40Ea  ;  Micro*  strissis:  sti  *  a,sto 
5841  PRINT  USING  a37X,7A,9.4DEa  ;  astl»  *  \stipy 
S85I  PRINT  USING  a38X,22Aa  ;  a(*f< qiodratic  R  4  R*>>* 

5841  PRINT  I  PRINT  ‘ItWIWBIMWWMWMmWttNWttlWttlWWt 

5871  BISP  aUsc  (ROLL)  key  ti  fin  CRT  oityit.a„ 'Press  CONTINUE* 

5881  PAUSE 

5891  CLEAR  I  PAGESIZE  14  I  DISP  8  DISP  |  BISP 

59(1  BISP  aBi  yi«  wish  ti  cinpetc  a  nen  set  if  nincnts,  Ni,H2,H4,  index  t," 
5911  BISP  "and  irientatiin,  THETAtT  (NO  will  terninate  prtqran  iperatim).' 
5921  BISP  *  Y/N" 

5931  INPUT  Al 

5941  IF  UPC«  (A*)«aNa  THEN  GOTO  4111 

5951  PRINT  I  PRINT  t  PRINT  *Inp»t  new  liads  Hl,N2,H4a  •  PRINT  I  PRINT 
5941  INPUT  Hi,N2,N4 

5971  PRINT  •  PRINT  aNw  nments:  N1  =a;«lja  H2  =*}«}•  N4  *a;H4 
598B  WAIT  INI 
5991  GOTO  5181 

4188  PRINT  •  PRINT  'End  tf  priqro*.  Thank  yie.* 

4118  OS 
4(28  ! 

4838  ! 

4(48  THETA:  !  Sukriotine  t»  enter  THETA'S  fir  cilcilatiin  if  Vi's  and  ViB's. 
4(58  FOR  t»ti  TO  n/2+ti-i 

4848  PRINT  8  PRINT  *  THETAta;t;a)  *  ?a  8  PRINT 

4(78  INPUT  THETA(t) 

4888  PRINT  ,THETA(t) 

4898  THETA(t)>THETA(t)«Pl  /1I8 

41(8  GOSUB  Vi 

4118  GOSUB  ViP 

4128  NEXT  t 

4138  x*28ho 

4148  Vl*xWl  8  V2*xW2  8  V3«xW3  8  V4»xlV4 
4158  i*2/3>i)ik3 

4148  UlBntVlB  8  V2>t«V2B  8  V3B*xtV3>  I  V4B*xlV4B 
6178  RETURN 


6181  ! 

6191  Vi: 

6211  Vl*Vi+CQS  (2»THETA(t» 

6211  V2«V2tC8S  (4*MTA(t)) 

6221  V3*V3+SIN  (2»THE7A(t)) 

6231  V4»V4*SIN  (4ITHETA(t)) 

6241  RETURN 
62SI  ! 

6261  ViD: 

6271  i«t*3-(t-i)*l 

6281  V1B*V1N*C8S  <2»THETA(t))«i 

6291  V2D*V29HSS  (4«THETA<« )  >lx 

6381  V3J-V3HSIN  (2>THETA( T) >8x 

6311  V4W4D+SDI  (48THETA(t))tx 

6321  RETURN 
6331  ! 

6341  ! 

6351  HUM:  !  TrvnftrMtiin  if  qindratic  itrcaqth  pvMttrs  in  strain  space 
6361  !  in  mltiplt  wqlt  fmctim.  17.6a,  p.298 
6371  mrfOS  (THETAt)  t  na«SIN  (TNETAt) 

6381  Cll<XX«MiA4iGYYtna*44€XYf2tM*2»M<24CSSt4tMi42lnn>2 

6391  C22^XX8HnA4*6n<mA4iGXYt2tm‘2tn«424fiSSt4tM<2lnn42 

6411  C12>CWt«i',2ton*2+€m«n‘2tefl*2+CX1f*(«i*4Hin*4)-<SS44»*m*2tnn‘2 

6411  G66<XX*Mi*2fM*2+6YYtM*2tNi‘2-CXYt2tM*2tAn*2*C$St(M*2-nn*2)*2 

6421  C16aQ(X>Mi<38im’6YYtMitnn*3'iCXYt(iMlM*3~M*3tM ) +CSSt2S(  iwtnn  *3-w»*3*nn ) 

6431  G26»€XXIni»lnn*3*CYYtitnA3lnn+0(Yt<nii*3lon-M*nn43)*5SSS<ct(nii‘3tnn-niilnn*3) 

6441  Cl>CXtiM*24€YfMi‘2  !  T7.66,  p .291 

6458  G2HX$n»*2+GYtimA2 

6468  C6~tXliin8wn  CYtnntni 

6478  a4ilttli*2+C22ta2i*2+28Gi2ttiit«2i+G66lt6i*2+2RCi6<Cli8t6i+28S26S«2tte6i 
6481  b*C18ali+C28t2i*G68c6i 
6491  c*”l 

6588  Rn*(-b+S8R  <b*2-«afc))/(2*n) 

6518  Rqp*AIS  ((-b-SOR  <b*2-4talc))/(2ta)> 

6528  RETURN 
6538  ! 

6548  ! 

6558  SSSR:  !  Subwtint  ti  cmpitt  mi  stress-strain  R  and  R* 

6568  Res*S/AIS  (esilES) 

6578  IF  (U>8  THEN  Rei»X/<exilEI>  ELSE  Rei«XP/A8$  <t*i«EX) 

6588  IF  ifi)8  THEN  R«f*Y/(«y*8CY)  ELSE  R«f*YP/AiS  UfotEYI 

6598  IF  RcxUev  THEN  Rt>Rti  ELSE  Re*Rej 

6688  IF  R«(Rt  THEN  Re«Res 

6618  IF  «u<8  THEN  lti*X/ABS  (eulEX)  ELSE  Rei*XP/(eii*EX) 

6628  IF  e»<8  THEN  R«f«Y/AIS  (ejilEY)  ELSE  Re»*YP/(t»*IEYI 
6638  IF  Rci<Rtf  THEN  Rep*«  aSE  RtHif 
6648  IF  Res(Rtp  THEN  Rcp*8et 

6658  Rss*S/AIS  (ssi) 

6668  IF  su>8  THEN  8si'X/su  aSE  Rsx*ff/ABS  (hi) 

6678  IF  Sfi)8  THEN  Rsf*Y/spi  aSE  Rsp«ff/AK  (svit 
6688  IF  RsKRsp  THEN  ls«Rsx  aSE  Rs«ts? 
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